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In this paper, we consider a Unear quantum network composed of two distantly separated cavities 
that are connected via a one-way optical field. When one of the cavity is damped and the other 
is undamped, the overall cavity state obtains a large amount of entanglement in its quadratures. 
This entanglement however immediately decays and vanishes in a finite time. That is, entanglement 
sudden-death occurs. We show that the direct measurement feedback method proposed by Wiseman 
can avoid this entanglement sudden-death, and further, enhance the entanglement. It is also shown 
that the entangled state under feedback control is robust against signal loss in a realistic detector, 
indicating the reliability of the proposed direct feedback method in practical situations. 

PACS numbers: 02.30.Yy,03.67.Bg,42.50.Dv 



I. INTRODUCTION 

Reliable generation and distribution of entanglement 
in a quantum network is a central subject in quantum 
information technology especially in quantum com- 
munication 0, IE 14, [i| ■ The biggest issue in such sys- 
tems is the decay of entanglement due to decoherence ef- 
fects that are inevitably introduced when node-channel or 
channel-environment interaction occurs. Entanglement 
distillation [1, 0| is a useful technique that restores such 
degraded entanglement. However, it sometimes happens 
that entanglement completely disappears in a finite time, 
which is called entanglement sudden- death [1, [^. In this 
case, distillation techniques cannot recover the vanished 
entanglement. 

On the other hand, feedback control can be used to 
modify the dynamical structure of a system and improve 
its performance, e.g., see [13, [ll!, [ll, [ll]- Entangle- 
ment protection or generation is one of the most attrac- 
tive applications of feedback [3, [H, [l^, In par- 
ticular, two studies have demonstrated that a feedback 
controller effectively assists in the distribution of entan- 
glement in a quantum network. One such result is by 
Mancini and Wisema n llSl . where a direct measurement 
feedback method [l^, l2d| is used to enhance the corre- 
lation of two bosonic modes that couple through a 
nonlinearity. The other such result is by Yanagisawa [2l| , 
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where an estimation-based feedback controller is used to 
deterministically generate an entangled photon number 
state of two distantly separated cavities. 

This paper follows a similar direction to [l^ and [2l| . 
That is, we also consider a problem of distributing en- 
tanglement in a quantum network using direct feedback 
control. The quantum network being considered is de- 
picted in Fig. 1: Two spatially separated cavities are 
connected via a one-way optical field, and the measure- 
ment results of the output of Cavity 2 are directly fed 
back to control both cavities. A more specific descrip- 
tion will be given in Sections II-B and ITC, but here we 
note that the network model considered brings together 
the following three features that have been not simulta- 
neously considered in previous work. First, the network 
contains models of realistic components; a realistic quan- 
tum channel in contact with an environment and a realis- 
tic homodyne detector with finite bandwidth ,2^, i2^, 2^ ■ 
A realistic model is of practical importance for real-time 
quantum feedback control. Second, we consider linear 
continuous- variable cavity models (i.e., we consider the 
quadratures of the cavity mode), similar to the case of 
0, H, [ill- Hence, the system differs from a discrete- 
variable system such as an atomic energy level [1, [3] , or 
a photon number system description [21]. This setup is 
motivated by the recent rapid progress and deep under- 
standing of continuous-variable systems in the quantum 
information regime [25^. Third, the cavities are spatially 
separated, and the interaction between them is simply 
mediated by an optical field, while in [ll| two bosonic 
modes interact through a x*-^^ optical nonlinearity and 
thus the two modes are not spatially separated. The spa- 
tially separated case is the case of interest in applications 
such as quantum communication. 
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FIG. 1: Schematic of the network. Thick gray lines repre- 
sent quantum optical fields, while thin black lines represent 
classical signals. 



The contributions of this paper are as follows. First, 
we show that the network considered in this paper, which 
looks complicated, can be systematically captured and 
described using the theory of quantum cascade systems 
[2l,[23,[2i,[2i]. We then show that when the first cavity 
is undamped and the second cavity is damped, the cav- 
ity state obtains a large amount of entanglement, which, 
however, disappears in a finite time despite the contin- 
uous interaction between the cavities; i.e., entanglement 
sudden-death occurs. As mentioned above, no distillation 
technique can recover such a zero entanglement. Never- 
theless, we show that direct feedback control not only 
prevents entanglement sudden-death, but can also en- 
hance the entanglement. Moreover, it will be shown that 
the entangled state under control is robust against signal 
loss in a realistic detector, implying the reliability of the 
direct feedback method in practical situations. 

We use the following notation. For a matrix A — {uij), 
the symbols , A^ , and A* represent its transpose, con- 
jugate transpose, and elementwise complex conjugate of 
A, i.e., AT = (a^,), At = (a* J, and A* = [a*^) = {A^Y , 
respectively. These rules are applied to any rectangular 
matrix including column and row vectors. Re(A) and 
Ini(A) denote the real and imaginary part of A, respec- 



tively, i.e., {R,e{A))i 



a*,)/2 and (Im(A)), 



[ttij — a*j)/2i. The matrix element can be an opera- 
tor hij] in this case, a*^ denotes its adjoint operator. 



II. MODEL 



A. General linear quantum systems 
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Suppose that the system contacts with M optical vacuum 
fields without scattering. In general, such an interaction 
is described by a unitary operator that obeys the Hudson- 
Parthasarathy equation [30|: 



dUt = 



M 



k=l 



M 

E 

fe=i 



LkdBlt ~ LldBk,t) Ut, Uo = i- (1) 



The operators Bk^t and B^ ^ represent the quantum anni- 
hilation and creation processes on the k-th field, respec- 
tively. Note that [dBi,t,dB* ^] = Sijdt. Let us choose 

H = H* and Lk as follows: 

iJ-ii^gi, Lk^Llk, (2) 



where G = G" € R^^vx^iv ^^^^ ^ ^zn _ rpj^g 
system variables obey the Heisenberg equation Xt := 
(. . . , iltqiUt, U^PiUt, ■ ■ y. We then obtain the follow- 
ing linear equation: 

dxt ^ Axtdt + lY.NiL'^dB; - L^Bt], (3) 

where {Li,...,Lm) e C^^x*^ A := Ew[G + 

lm{L^L)], and Bt := {Bi^t, ■ ■ ■ , BM,tV ■ For details 
on the physical meaning of these abstract linear mod- 
els, see Section III and, e.g., [sil, [H, [s^l. It is 
easy to see that the first moment vector (xt) := 
(. . . , {U;q,Ut), {Ut*P^Ut), . . y, where (x) := Tr (pi), sat- 
isfies the linear equation d(xt)/dt = A{iit). Also, the 

covariance matrix Vt = ( ( Vy ) ) , where 



Ax(AxJ -I- (AxjAx 



Ax, 



satisfies the following Lyapunov matrix differential equa 
tion: 

dVt^ 
dt 



= AVt + VtA^ D. 



(4) 



Here, D := SjvRe(LtL)S]J^. Suppose that the quantum 
state is Gaussian at i = 0. Then, from the linearity of 
the dynamics, the unconditional state is always Gaus- 
sian with mean (xt) and covariance Vt. Note that the 
unconditional state corresponds to a classical probability 
density that describes a linear diffusion process. 



We consider a general linear continuous-variable sys- 
tem with A^-degrees of freedoms. Let x^ = {qi,py be 
the standard quadrature pair of the i-th subsystem. It 
then follows from the canonical commutation relation 
[SiiPj] = QiPj ~ PjQi = i^ij = 1) that the vector of 



B. The ideal network 

Before describing the quantum network depicted in 
Fig. 1, let us consider the ideal situation where the opti- 
cal field between the cavities is not in contact with any 
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environment, and the homodyne detector is perfect. In 
this case, the system is the simple cascade of two cavities 
with a feedback loop. The entangled state of this ideal 
network will be compared to that of the realistic one for 
the purpose of clarifying how much the realistic parame- 
ters affect the system. Also, this ideal setup allows us to 
determine a reasonable control Hamiltonian (the vector 
/ given below), as will be seen in Section III-B. 

Each component of the network is described as fol- 
lows. The optical vacuum field Bi^t comes into Cavity 1, 
and then, its output becomes the input of Cavity 2. We 
assume that, after some approximations, the i-th cavity- 
field interaction is represented by Eq. ^ with single field 
(M = 1) and with the following operators: 



where = (di-|-a*)/\/2 andpi = (di — a*)/\/2i. The sub- 
script (1, i) in L means the 1-st field and the i-th cavity 
(see the figures in Appendix A). Also, d ^ Gj € M^^^ 
and £i e C^. The output of Cavity 2 is transformed to 
a classical signal yt through an ideal homodyne detector. 
Suppose now that each cavity has an additional Hamil- 
tonian of the form 



fb 



utFi^utfJ^i, fie 



1,2) 



where G R is the control input. Then, direct measure- 
ment feedback ut = gyt closes the loop by connecting 
the detector to the cavities. Here g e M is the control 
gain. Note that we need a classical communication chan- 
nel in order to control Cavity 1; that is, a local operation 
via classical communication (LOCC) type control is per- 
formed. 

For this network, we can easily determine the system 
matrices G and Lk in Eq. ([2]) that specify the whole 
dynamical equation. The derivation is based on the the- 
ory of quantum cascade systems [1^, [13, [H, [1^ and is 
given in Appendix A. Then, from the definition, the A 
and D matrices in the Lyapunov equation (|4]) are readily 
obtained as follows: 

Aid = Ao + 2gE2/Re(£)T, (5) 

Ad = Do + Md^ff^ - 9iMe)f - gfiMeV]^l, (6) 

where £ = (4, 4)^, / = ifl , fl)\ and 



Ai 
2SIm(^;4) A2 



Di * 
ERe(£^^]")ET D2 



(7) 

Here, A, = E[G, + lYn{£*£j)], Di = Y.Re{£* £j , and 
denotes the symmetric elements. Note that Aq and Do 
are the system matrices of the network without feedback. 
Hence, the upper off-diagonal block matrix of Ao is zero, 
implying the one-way interaction of the cavities. 



C. The realistic network 

We are now in the position to describe a realistic net- 
work, which introduces the following two assumptions. 
First, the output of Cavity 1 is mixed with another vac- 
uum field ^2,4 through a beam splitter (BS) with trans- 
mittance a. This is a standard model of possible envi- 
ronmental effects on a long quantum channel. Second, 
the homodyne detector is not perfect and is described by 
the one-dimensional classical dynamics 

d^t = aiitdt + a2dwt, dyt = az^tdt + dvt, a, e K, (8) 

where wt is an input stochastic process satisfying 
E[du;(] = dt and vt is an additional measurement noise 
satisfying E[c?W(] — a^dt (04 > 0). In particular, atypical 
low-pass filter (LPF) is realized by choosing as 



ai 



1 1 
T r 



as 



where r > is the time-constant. We now note that 
the detector ([5]) can be represented as a quantum system 
with two fields wt — Bn + and vt = Bs^t + ^3 t, 
where Bi ^ is the output of Cavity 2. Indeed, from [29j . 
Eq. ^ with M ~2 and with the operators 



H3 = y(g3P3 



P393), ^1,3 = -102^3, ^3.3^7, 93 

2a4 



leads to a linear equation of the form ([5]), where Ufq3Ut 
plays the same role of . 

Consequently, the network is composed of two cavities, 
a beam splitter, a detector, and a controller, with three 
optical vacuum fields. (Note that the beam splitter with 
local oscillator (LO) shown in Fig. 1 is a part of the 
detector.) To systematically obtain the overall system 
matrices G and Lk in Eq. ([2]) of this complicated network, 
we again use the theory of quantum cascade systems [1^, 
12^ . [23, [2^. The procedure is given in Appendix A. We 
then obtain the matrices A and D in Eq. (|4|) as follows: 



Ai-c — 



Ai 

2aY:im{£*2£j 



503 S/i 
A2 503^/2 
2aa2Re(£i)"r 2a2Re(^2)^ ai 

Di ★ 
aSRe(^^4)ST D2 
-aa2lm(£i)TET -a2lm(^2) 

+ <?^a4 S/2/7ST ^f2fj^'^ 




(9) 



^TvT 



(10) 



III. ENTANGLEMENT CONTROL 

In this section, we study the entanglement of the cav- 
ity state for the ideal network. Since the state is Gaus- 
sian, its entanglement is completely characterized by the 
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covariance matrix [3J, |35|. In our case, the covariance 
matrix to be evaluated is obtained by solving the Lya- 
punov equation (j4|) with the coefhcient matrices Aid and 
Did in Eqs. ^ and ([5]). Let the matrices Vi be defined 
by the 2x2 block matrix decomposition of V as follows: 



V = 



Vl V2 



Then, the following logarithmic negativity |36| can be 
used as a reasonable measure of Gaussian entanglement: 

^AA^maxjO, -log(2z/)}, (11) 

where logo; denotes the natural logarithm of x, 



i.:=-^yA-.^A2-4dct(l/), 
A := det(V"i) + dct(y3) - 2det(V2). 



(12) 



The logarithmic negativity £V divides the state space 
into two regions: (i) the separable region, corresponding 
to Ej^ = 0, and (ii) the entangled region, within which 
£V > 0. Thus phenomena of entanglement creation and 
destruction can be understood simply in terms of move- 
ment of the system between these two regions. 



A. Entanglement sudden-death 

Here we study the uncontrolled network; i.e., g = 0. 

We compute Ej\/ for two situations. First, consider the 
case where both cavities have the same quadratic Hamil- 
tonian and are damped as a result of the field-cavity in- 
teraction, that is. 



Gi — G2 



m 
1 



where m > and k > 0. This type of quadratic Hamilto- 
nian can be implemented in a cavity system followin g th e 
scheme of the degenerate parametric amplification [28| : 
see Appendix B. In this case, Ai^ is a stable matrix, and 
the Lyapunov equation ^ has a unique steady state so- 
lution (see e.g., (s^)- Now assume that at t = 0, the 
cavity is in the separable state satisfying Vq = 21. When 
the optical field is switched on, the cavity modes couple 
after a finite time (i.e., the "entanglement sudden-birth" 
[sst occurs) , and a steady entangled state is generated as 
seen from the dotted line in Fig. 2 (a). However, in this 
case, the entanglement is very small {Ej^ w 0.21). 

This result can be understood by examining the trajec- 
tory of the parameter (A,det(V^)). In Fig. 3, the colored 
region with contour lines represents the set of parameters 
where a general two-mode Gaussian state is entangled, 
i.e., Ej^ > 0, while the white region corresponds to sep- 
arable states; i.e., E_\f = 0. The trajectory, denoted by 
Tdaim evolves toward the steady entangled state that is 



located far from the area with large Ej^ (the right bot- 
tom area in Fig. 3). This is likely because each cavity has 
a strong tendency to transit into the vacuum state due to 
the damping. Indeed, when the cavity is in the separable 
vacuum state |0) |0), the corresponding covariance matrix 
satisfies (A, det(F)) — (0.5,0.0625), which is very close 
to the equilibrium point of Tdam- Moreover, Fig. 2 (b) 
shows that the purity (for a discussion of physical mean- 
ings of the purity, see e.g. [3§|) of the steady Gaussian 
state p, 



Tr(/52 



1 



e (0,1], 



(13) 



approaches PwO.Sasf— >cx). This also suggests that 
the steady state is close to the separable vacuum state. 

The above observation motivates us to try a dispersive 
field-cavity interaction, which results in a phase shift of 
the output field [13, Eo, SH IH ■ For a practical method 
to implement this kind of coupling in a cavity system, 
see Appendix B. In this case, the cavity is not damped, 
and thus, it does not have a tendency to move toward 
the vacuum state. In particular, we assume that only 
the first cavity has such an interaction; i.e., 

h = iV^, 0)T, £2 = y^(l, i)^. 

We then find that A-^i is not a stable matrix, and the Lya- 
punov equation ^ need not have a steady state solution 
as t — > 00 [131 • Fig. 3 shows that the corresponding tra- 
jectory, denoted by Tdis, evolves far from the separable 
initial state and reaches the area with large Ej^. This 
figure also shows how both the entanglement and purity 
decreases as time goes on and Tdis escapes from the region 
of entangled states at t = 6.2. 

Finally, we remark that, if we exchange the order of 
the interactions, i.e., £1 — y/K{l, i)^ and £2 = (v^, 0)^, 
the corresponding trajectory remains within the region of 
separable states, i.e., Ej^ = for all t > 0. The situation 
is much the same even when each cavity interacts with 
the field in a dispersive way. 



B. Feedback control 

We first discuss how to determine the coefficient vec- 
tor / — {fj, fJY that realizes high-quality entanglement 
control. Fortunately, in the ideal case, we can explicitly 
find such an /. The idea was originally provided by Wise- 
man and Doherty in [43| , but here we apply the idea in 
a slightly different manner. 

Assume that g = 1. Then, the Lyapunov equation ^ 
with coefficient matrices Aid and I?id in Eqs. ([5]) and ^ 
can be rewritten as 



dt 



^niVt) + Mf-ftKf-ftV^l, (14) 



where 



ft := 2J:2VtRe{£) + lm{£) 



(15) 
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(a) 



(b) 



FIG. 2: Time-dependence of the logarithmic negativity (a) 
and the purity (b) of the overall cavity state without feed- 
back control. The solid and dotted lines correspond to the 
dispersive-damped and damped-damped cases, respectively. 
The parameters are m — 0.2 and k = 1. 




FIG. 3: Trajectories of the parameter {A,det{V)) without 
feedback control. Tdam and Tdis correspond to the damped- 
damped and dispersive-damped cases, respectively. We set 
Vo = 2I a.tt = 0, from which A = 8 and det{Vo) = 16 follow. 
(A and det(l/) are in units of ft = 1 and — 1, respectively.) 



and 

niV) := AoV + VAl + Do 

2VRc{e) - S2lm(£)j \2VRe{£) ~ E2lm(^) 



We now recall from Fig. 2 (b) that entanglement sudden- 
death occurs simultaneously with a decrease of the purity 
(fTS]) . This suggests that preventing a decrease of purity 
may also prevent entanglement sudden-death. However, 
we should point out that it is not apparent that this will 
always be the case and the relationship between loss of 
purity and entanglement sudden-death needs to be stud- 
ied further. Therefore, a simple control strategy that we 
try here is to find a feedback controller that prevents an 
increase of det(Vt) in order to keep the purity high. As 
the second term of the right-hand side of Eq. is 
always non-negative, it is then reasonable to choose the 



time- variant coefficient vector f = ft- Of course, this 
intuitive argument does not always allow us to conclude 
that det(Vt) takes its minimum value. However, it is 
known that the algebraic Riccati equation TZ(V) ~ has 
a solution satisfying det{V) — 1/16, which implies that 
the maximum purity P = 1 is achieved; e.g., see [43l |. 
Now assume that Eq. (fTi|) has a unique steady solution 
Voo for a constant /. Then, by taking the time- invariant 
coefficient vector 



/ := 2S2l4oReW + Im(£), 7^(Foo) - 0, 



(16) 



we obtain the same desirable result, det(V|3o) = 1/16. 
Note that the numerical solution to the algebraic Riccati 
equation 7?,(V"oo) = can readily be obtained using a 
standard software package such as matlab. 

We now consider direct feedback control with the co- 
efficient vector ((T6)) . Let us begin with the case where 
the first cavity-field interaction occurs dispersively. For 
this system, it is expected from Fig. 3 that the trajec- 
tory Tdis can be modified and stabilized via feedback so 
that it has an equilibrium point in the area where Ej^ 
is large. That this is indeed true is shown below. When 
the parameters are given by m = 0.2 and k — I, we find 
that / = (0.1212, 2.2196, -0.3163, -3. 2277)T from (fT6|). 
Fig. 5 illustrates that the controlled trajectory, denoted 
by 72is, indeed shows the convergence that we had hoped 
for. The entanglement and the purity of the steady cav- 
ity state are shown in Fig. 4. While / is determined with 
fixed g = 1, we consider variations in g to gain under- 
standing about its effect on the control system. When 
control is not used {g — 0), the pair of dispersive and 
damped cavities does not settle down to a steady state 
as seen in Section HTA, and we find that E_\f ^ as 
i — > cxD. On the other hand, even with the small-gain 
feedback controller, the system becomes stable and has 
a unique steady state with nonzero entanglement. Re- 
markably, when (7 = 1, the entanglement of the steady 
state {Ej^ w 2.2) improves upon the maximum value 
of E_\f of the uncontrolled state {Ej^ « 0.65) shown in 
Fig. 2 (a) . Hence we see that direct feedback not only pre- 
vents entanglement sudden-death, but can also enhance 
the entanglement. 

Feedback can also improve the entanglement of a sys- 
tem where both cavities are damped, but it is still very 
small as seen from the dotted line in Fig. 4 (a). (The 
coefficient vector defined by Eq. (fTB]) in this case is cal- 
culated to be /= (0.0629,0.1525,0.2479,-0.5830).) To 
understand this phenomenon, we recall that the uncon- 
trolled trajectory Tdam has an equilibrium point that is 
located far from the area with large Ej^. Hence, it should 
be hard to drastically modify this trajectory such that it 
could reach that area. It is actually observed in Fig. 5 
that, even with the vector /, the controlled trajectory 
'^'am shows almost the same time-evolution as the au- 
tonomous one Tdam- 

The above results suggest that strong stability of the 
autonomous system sometimes makes it difficult for the 
state to transit into a desirable entangled target. 
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FIG. 4: The logarithmic negativity (a) and the purity (b) 
of the steady cavity state with feedback control, g is the 
control gain. The solid and dotted lines correspond to the 
dispersive-damped and damped-damped cases, respectively. 
The parameters are m = 0.2 and «: = 1. 
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FIG. 5: Trajectories of the parameters (A,det(V)) with 
feedback control. T^i^^ and T^^g correspond to the damped- 
damped and dispersive-damped cases, respectively. The ini- 
tial state is the same as before: Vb = 27. 



IV. A REALISTIC CONTROL SCENARIO 

Finally, we return to the original setup of the network. 
That is, the quantum channel is in contact with an en- 
vironment, and the homodyne detector is replaced by a 
realistic LPF with finite bandwidth. The purpose here is 
to study the impacts of these realistic components on 
the entanglement of the cavity state. The covariance 
matrix of the cavity state corresponds to the left-upper 
4x4 submatrix of Vt that is the solution of Eq. ^ with 
Aj-c and D^c given in Eqs. ([9]) and pO)l . Note that the 
cavity state is the reduced one with the detector mode 
traced out. We here focus only on the network where 
the first cavity interacts with the field dispersively. For 
the control Hamiltonian, we use the same coefficient vec- 
tor / = (0.1212, 2.2196, -0.3163, -3.2277)^. It should be 
noted that, in this realistic case, we cannot follow the dis- 
cussion in Section III-B to obtain a reasonable coefhcient 



FIG. 6: The logarithmic negativity of the steady cavity state 
with feedback control, g is the control gain, (a) From the 
top downwards, the lines correspond to r = 0.01, 0.2, 0.4, 0.6, 
while a — 1. (b) From the top downwards, the lines corre- 
spond to a = 1, 0.99, 0.95, 0.90, while r = 0.01. In both cases, 
the LPF noise is very small; 04 = 0.01. The parameters are 
m = 0.2 and k — 1. 



vector /. 

First, consider Fig. 6 (a). This shows some plots of 
Ej^ with the time-constant r changing between 0.01 < 
r < 0.6 and with the fixed transmittance a = 1 (i.e., no 
loss in the channel). The most upper line almost coin- 
cides with the ideal one shown in Fig. 4 (a). That is, the 
entanglement in the realistic situation continuously con- 
verges to the ideal one as r ^ 0. We also observe that 
the degradation of is small with respect to t. Since 
the detector is regarded as a component of the controller, 
these results imply that the realistic direct feedback is ro- 
bust against signal loss in the LPF. In other words, direct 
feedback control is reliable even in this realistic situation. 

On the other hand, Fig. 6 (b) plots Ej^ for some values 
of the channel loss f3 := 1 — a with fixed r — 0.01. We 
find that Ej^ converges to the ideal one as /3 ^ 0, simi- 
lar to the above case. However, in this case, E^ rapidly 
decreases with respect to (3. Even for the very small loss 
f3 = 0.01, a visible degradation occurs. Moreover, when 
f3 = 0.1, which still means we have a high-quality quan- 
tum channel, Ej^ decreases less than half of the ideal one. 
That is, the entanglement is fragile to realistic channel 
loss. 

The above results are reasonable because the channel 
loss directly reflects the decrease of interaction strength, 
while the finite bandwidth of LPF simply implies loss of 
a classical signal. Hence the former should be a critical 
factor for entanglement generation. 



V. CONCLUSION 

The contributions of this paper are summarized as fol- 
lows. First, it was shown that, when the first cavity is un- 
damped and the second one is damped, the overall cavity 
state obtains a significant amount of entanglement, which 
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however disappears in a finite time. Then, we have shown 
that direct measurement feedback can avoid this entan- 
glement sudden-death, and further, enhances the entan- 
glement. Moreover, it was shown that the direct feedback 
controller is reliable under the influence of signal loss in 
a realistic detector, although imperfection in the quan- 
tum channel is a critical issue that largely degrades the 
achieved entanglement. We believe that the case study 
we have presented provides useful insights that may be 
of use for more complex quantum network engineering. 

APPENDIX A: QUANTUM CASCADE SYSTEMS 
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• 

1/2,1 
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• 

L2,2 






• 


B'M,t 


• 

Lm,2 





FIG. 7: Abstract illustration of the cascade system. The 
black circles represent that the subsystem interacts with the 
field. 



1. General theory 

In this appendix, we begin with a review of the theory 
of quantum cascade swtems which was originally devel- 
oped by Carmichael [2^, H^l and Gardiner [28j] in a quan- 
tum optics framework and recently reformulated in more 
general setting by Gough and James [1^ . We then apply 
the theory to our model and derive the corresponding 
system matrices. 

The most general form of quantum dynamics that in- 
teracts with AI optical input fields is described by the 
following unitary evolution: 



dUt 



M 



k=l 



M M 

'^LkdBl t- L*SjkdBj,t 

k=l k,j=l 
M 

{Sk] ~ 5k])dkkj 



Uu Uo^I. (Al) 



This is also called the Hudson-Parthasarathy equation 
[30| . Here, the operators Bk,t and Bl ^ represent the 
quantum annihilation and creation processes on the k-th 
field, respectively. The operator Akj represents the scat- 
tering process from the k-th state to the j-th state, and 
it satisfies dAijdAki = SjkdAu. The matrix of operators 
S := (Sij) must satisfy S^S = SS^ = / in order for Ut 
to be unitary. The system is completely characterized 
by the triple T = {S,L,H), where L is a vector of op- 
erators L := (ii, . . . ,L]\fy. Let X be an operator of 
the system. Then, this evolves in time according to the 
Heisenberg equation X jtiX) := tj^XUt- In particu- 
lar, we can define M output fields i?^ ^ :— jt{Bk,t), which 
yields 

dB[^jt{t)dt + jt{S)dBt, 

where we have defined Bj :— {B[ j, . . . ,B'j^^.nY,jt(L) '■= 
(jt(Li), ... ,jt(LM))"^, etc. 

Let us now consider two systems: Fi = (Si,Li,//i) 
and r2 = (S2, L2, -??2)- Note that the number of inputs 



(i.e., outputs) of these systems can always be matched 
by introducing additional components in L and J in S 
as L © and S /. These systems can be connected so 
that the outputs of Fi are the inputs of F2, as depicted 
abstractly in Fig. 7. We denote this cascade system by 
F2 < Fi. Then, from [H, we have 

F2 <Fi = fs2Si, L2 + S2L1, 



^i+^2 + ^(L^S2Li-L1s^L2)), (A2) 



where ^ = (LJ_„ . . . , Ll,^,), (k = 1, 2) and St = (S*,). 

Direct measurement feedback [l^ [23| is no more than 
a cascade of the system and the controller. Hence, the 
overall system representation of the controlled network is 
readily derived using Eq. (jA2[) as follows. For simplicity, 
let us consider a single-input single-output system F = 
{S, L^H -\- UfF), where Uf represents the control input. 
An ideal homodyne detector yields a classical signal yt = 
jt{Bt + Bl). Then, the direct feedback ut = gyt closes 
the loop and realizes 

Ffb = (l,-igF,0)<(5,L,i7) 

= (5, L - igP, H + ^{FL + L*F)'^ . (A3) 

For a more detailed discussion, see 12911 . 



2. Ideal network 

Let us then apply the above formulas to our system. 
First, we consider an ideal network composed of the fol- 
lowing three subsystems: 

Cavity 1: Ci - (/, Li,i, H^), 
Cavity 2: C2 = (/, Li,2, ^2), 
Controller: F = (/, -igP, 0), 

where Li^i and Hi are given in Section II, and F := 
Fi + F2 = fjxi + /Jx2 . The abstract configuration of 
the network is given in Fig. 8. Thus, the network is given 
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FIG. 8: Abstract illustration of the ideal network. 
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FIG. 9: Abstract illustration of the realistic network. 



by 



where 



and 



F < C2 < Ci = (/, Lid, i^id), 



Lid = + ^1,2 - ^gF 

= [£j - igfj, 4 - igfj] 



Xl 
X2 



^^id — Hi + H2 + ^(-^1,2^1,1 ^ Ll 1L12) 



+ 



lm{£i£. 



G2 



.g/Rc(£)T + 5Rc(£)/T}x 



=: iic'^GidX. 



(A4) 



F(il,l+il,2) + (i^ l+^t2)^^ 



Gi Im(£i4) 



Here, i := (4,4)"^ and / {fJjJV- From the def- 
inition, we then obtain the system A and Z? matrices: 
Aid := S2[Gid+Im(L|'^Lid)] and Ad := S2Re(Lj^Lid)ST, 
and they are given in Eqs. (O and 



3. Realistic network 

We next consider a realistic model of the network. 
Each component is given as follows. 



Cavity 1: Ci = /, 



Beam Splitter: B 








, Hi 



a -/? 
f3 a 
1 



0, 0), 



Detector: D = / 



^1,2 

Cavity 2: C2 = (/, , if^ 



, H. 




-igP 

The /c-th element of the above vectors corresponds to the 
k-th quantum field Bk,t- Fig. 9 abstractly illustrates the 
structure of the interactions in the network. Note that 
the detector part includes the beam splitter with local 
oscillator shown in Fig. 1 . Iteratively using Eq. (IA2[) , we 
obtain 



Controller: F = / 



F<]D<]C2<]B<]Ci 



a 




0" 




f3 


a 













1 





aLii + iyi_2 + -^1,3 
pLi 

L3.3 - igF 



H,, 



where 



Hrc — Hi + H2 + H3 + ^(-£1^2-^1,1 ~ -£1,1^1,2) 



1 

2i 



L 



{all , + it 2)-^!, 3 



3,3^)- 



It should be noted that 04 appears in the last term of H,-c 
due to dB3,tdB^^ = aidt. We now look at the relation 
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(S,L,iJ) = (S,0,0) < (/,StL,ff). This implies that 
any system (S, L, H) is equivalent to the system without 
scattering noises, (/, S^^L, H), as long as we focus on the 
unconditional state. This is because (S,0,0) just cor- 
responds to the modulation of the output that is to be 
discarded. (Note that if we consider the conditional state 
based on the measurement result, the above equivalence 
does not hold.) Consequently, our system is given by 



where 



Lii + ail, 2 + ail, 3 

^3,3 - igF 

Ij aij -iaa2 
-pe] i/3a2 

-^gfi -^gfl 03/204 



(A5) 









X2 




. ^3 . 



(A6) 



and Hr, 



Grc — 



x^G'rcx/2 with 



alm(£i4) 035/1/2 aa2Re{£i) 

G2 035/2/2 a2Re(^2) 
035/7/2 a35/2'"/2 ai 
aa2Re(£i)"r a2Re(£2)"^ oi 



(A7) 



Hence, we can now obtain the drift and diffusion matri- 
ces: 



Ac = S3[ac + Im(Lj,L,e)] = 

Ac - S3Rc(iJ,L,e)I]J - 



Arc 

-ai 



Drc 
03/404 



(A8) 



where A^c and D^c are given in Eqs. (O and PTI| . Since 
the 6-th column and low elements do not affect the others, 
it suffices to consider the reduced 5x5 matrices A^o and 

Drc. 



APPENDIX B: REALIZATIONS OF LINEAR 
STOCHASTIC MODELS 

The purpose of this section is to discuss possible phys- 
ical realizations of the linear stochastic models making 
up the entanglement scheme considered in Sections Iflll 
and If VI Linear models are approximations of real physi- 
cal system that are valid under various assumptions such 
as the dipole moment approximation and rotating wave 
approximation. In particular, we will describe approx- 
imate realizations of these models using optical cavi- 
ties. Discussions of the physical meaning of the abstract 
linear models considered herein can be found in, e.g., 
[3l|,[3l[^. 



1. Quadratic Hamiltonian 

Let a = {q + ip)/V2 and a* = {q — ip)/^/2 be the cavity 
annihilation and creation operators. Then, a quadratic 
Hamiltonian 



H = Aa*a + -{ee"^a*^ - e*e""^r) 



can be realized with a degenerate parametric amplifier 
(DPA) with a classical pump [2^, Section 10.2] in a ro- 
tating frame at half the pump frequency, where ee"^ (e, (f> 
real) is the effective pump intensity and A is the detuning 
frequency of the cavity mode of the DPA from the half 
the pump frequency (i.e., A = LUcav — t^p/2, where LOcav 
is the cavity resonance frequency and cUp is the pump 
frequency). It is then easy to verify that by choosing 



A = 



1 



1 



= 



the Hamiltonian can be written, in terms of the quadra- 
tures, a,s H = {m(p +p^)/2. The latter is the form of the 
Hamiltonian used in our models. 



2. Models with dissipative coupling L — ^/ko, and 
direct measurement feedback 



Linear systems with the dissipative coupling L = y/Jia 
are quite standard and can be implemented as an optical 
cavity with a leaky mirror, but here we shall also consider 
how the direct measurement feedback term F — utf^x., 
with / = (/i,/2)^, can be implemented in this system. 
Such an implementation is shown in Fig. 1101 The cav- 
ity has two partially transmitting mirrors Mi and M2 
with coupling constants ^/k and ^/j, respectively. Here 
7 is chosen such that 7 <C 1 and 7 ^ k. The cavity 
interacts with an incident vacuum noise field at mirror 
Ml via the dissipation coupling L — ^/na. The feed- 
back F is implemented as follows. First, the (real- valued) 
control signal Ut is amplified with gain 1 / ^/7 and mul- 
tiphed with / — ^(— /2,/i)^ to give the real signal 
ft = (I'l.t, i'2,t)''' = -^fut- vt is then sent to a modula- 
tor that displaces a vacuum bosonic field by the classical 
field vfds with = vi^t + i''^2.t to produce a coher- 
ent control field iic.t [Slj satisfying duc,t = v^dt + dBo,t, 
where -Bo,t is a vacuum noise field independent of -Bm.t- 
This displacement can be physically implemented by 
an electro-optical modulator, see |4^, Section III-B.5]. 
Mathematically, the displacement of a vacuum field Bg^t 
by a classical field w^ds is represented by the unitary 
Weyl operator M^(w^) (here wjj'(s) = for < s < i, 
and otherwise) satisfying the quantum stochastic dif- 
ferential equation (QSDE): 

1, 

t\ ■ 

F 



dW{vf{) = {vf[dB:^, - v?{dBo,t - '-\vf;\''dt)W{vf;)- 
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FIG. 10; Implementation of a dissipative coupling together 
with the linear control. 



with which we can write Uc,t = W{vf^)* Bo,tW{yf^). The 
coherent field Uc,t then interacts with the cavity via 
mirror M2 with coupHng coefficient ^^7, thus the to- 
tal cavity-fields interaction is described by the following 
QSDE: 



dUt 



"adt + \/7i{ddB*„ ^ — a*dBin,t) 



— — a*adt + y/^{ddu*^^ 



*dur 



For a sufficiently small value of 7, the effect of the noise 
Bo,t can be considered to be negligible and if also 7 <C k 
then its contribution to the system noise will be negligible 
compared to that of Bin,t- As a result, we find that the 
feedback term is included in the interaction: 



dUt = 



— iFdt — —a* adt + \fK{adB*.^ ^ 



a*dBr 



Ut 



The entire scheme is depicted in Fig. [TOl Note that a 
pumped x*-^^ nonlinear crystal can be placed inside the 
cavity to implement these linear couplings together with 
the quadratic Hamiltonian discussed in the preceding 
subsection. 



Models with dispersive coupling L — ^/kq and 
direct measurement feedback 



For realization of a dispersive coupling of the form 
L = \/Kq, consider the configuration shown in Fig. 1111 
This configuration consists of a ring cavity with mode a, 
an auxiliary ring cavity with mode b, a x^^^ nonlinear 
crystal in which the cavity modes a and b interact with a 
classical pump beam, and a beam splitter. The frequency 
of the auxiliary cavity is matched to half the frequency 
of the pump beam. The classically pumped nonlinear 
crystal implements a two-mode squeezing Hamiltonian 
given by Htms = ^(ee^''^p*a*S* - e*e'"p*a6) (where e is 
the effective intensity of the classical pump and LUp is the 
pump frequency), while the beam splitter implements the 
Hamiltonian Hbs = aa*b + a*ab* for a complex param- 
eter a. Suppose now that e and a are chosen to satisfy 
e/2 = ia = F for a real constant F > (in particular 
a — — iF). Then in a rotating frame at the frequency 
i-iJp/2, the overall interaction Hamiltonian between the 
modes a and b is thus given by 



i 




FIG. 11: Configuration of two cavities, a two mode squeezer 
(depicted by the square with an arrow to indicate classical 
pumping), and a beam splitter, to implement a dispersive 
coupling of the cavity mode a to the field Bin,t when the 
second cavity mode b is adiabatically eliminated. Note that 
Bin^t is phase shifted by 180° before reaching the mirror M. 



Assume that the coupling coefficient of the mirror M 
is large such that the mode b is heavily damped com- 
pared to mode d and has much faster dynamics than a. 
For simplicity, in the following we will use the formalism 
of quantum Langevin equations and a formal method to 
show that the configuration shown implements the dis- 
persive coupling L in the reduced dynamics for mode o 
only , after mode b is adiabatically eliminated [i^. See 
[33 | for a more rigorous derivation using QSDEs and the 
mathematical theory for adiabatic elimination developed 
in dl. 

Let Bin^t be an input field and Bout.t be an output 
field coupled to b at the mirror AI as shown in Fig. [11] 
and suppose that M has the coupling coefficient 7. In 
particular, notice the 180° phase shift in front of Bin,t 
before it strikes the mirror. Define r) to be a quantum 
white noise formally related to Bin,t as Bin^t = /q fi{s)ds 
and flout be the output noise at the mirror M formally 
related to Bout,t by Bout,t = Jq flout{s)ds. The quantum 
Langevin equations for the dynamics of a, 6, and f/out are 
[H, Chapter 5] 



b = 



F(6*-6), 

'^[HabM-^b + ^ri 



T{d* 



7 



Vout = y/^b-T]. 



(Bl) 

(B2) 
(B3) 



iT{a*b* — db — a*b + db*). 



Setting 6 = and solving Eq. (|B2p for b in terms of d 
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a and fj we obtain 

h = - {Via* + a) + V777) . (B4) 
7 

Substituting Eq. (|B4l) into Eqs. ((BT|) and (|B3)) we obtain 
that the reduced dynamics for a only and fjout are given 
by the following quantum Langevin equations: 

a = V{b*-h) 

2r 

= —{-V + r). (B5) 

2r 

?7out = — (a + a*)+77. (B6) 
The pair of equations (|B5p and (|B6p shows that the re- 



duced system after adiabatic elimination of the mode 6 is 
a single degree of freedom quantum harmonic oscillator 
with mode a coupled to the field i?m,t by the linear cou- 
pling operator L = 2y/2Tq/ where q — 
producing the output field _Bo„t^t = Jq fiout{s)ds. By suit- 
ably choosing T and 7 such that 2^/2T = yfRrf we see 
that with this scheme it is possible to implement any 
dispersive coupling of the form L = ^JHq. Moreover, 
by placing a pumped nonlinear crystal inside the cav- 
ity (pumped with the same frequency Wp) and adding a 
partially transmitting mirror in the ring cavity of a that 
couples it to a control field, one can easily combine this 
dispersive coupling together with the quadratic Hamil- 
tonian in Subsection 1 of this Appendix as well as the 
control shown in Fig. [TUl 
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